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Computer  Solutions  to  Heat  and  Diffraction 
Equations  in  High  Energy  Laser  Windows 

Volume  I 


1.  INTRODUCTION 


When  a  nonuniform  high  power  laser  beam  traverses  a  nearly  transparent 
material  medium,  the  imall  amount  of  energy  which  is  absorbed  by  the  solid  is 
often  sufficient  to  raise  its  temperature  by  an  appreciable  amount.  Since  this 
heating  is  also  nonuniform,  it  creates  thermal  stresses  throughout  the  window 
which  causes  distortion  and  defocusing  in  the  transmitted  light  beam;  this  phenom¬ 


enon  i3  called  thermal  lensing.  This  degradation  of  beam  energy  has  been  recog¬ 
nized  as  one  of  the  principal  failure  modes  in  high-power  laser  windows. 


In  order  to  assess  thermal  lensing  in  an  analytical  manner,  one  must  first 
determine  how  the  temperature  is  varying  throughout  the  medium  as  a  function  of 
time.  This  information  must  then  be  incorporated  into  a  formalism  which  charac¬ 
terizes  an  opti.cal  beam  traversing  a  material  medium  experiencing  stress 
birefringence  effects.  From  this,  it  now  becomes  possible  to  calculate  how  the 


Received  for  publication  26  November  1976 
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transmitted  beam  will  vary  with  time  in  the  vicinity  of  the  distant  Gaussian  focal 

.  ,  2-7 
point. 

To  quantify  the  effects  of  thermal  lensing,  an  efficient  computer  program 
package  has  been  developed  at  AFCRL  and  programmed  to  run  on  a  CDC6600 
computer.  The  package  has  been  written  to  handle  an  unpolarized,  symmetric, 
Gaussian-shaped  beam,  incident  uniaxially  on  either  a  thin  disc  —  or  annular¬ 
shaped  cylindrical  window,  which  is  assumed  to  be  homogeneous  and  perfectly 
antireflective  coated.  Three  coupled  programs  make  up  the  package.  One, 
called  TEMP5,  solves  the  full  heat  transport  equation  within  the  window  for  any 
given  set  of  initial  and  boundary  conditions  on  each  surface.  Its  output  is  fed  into 
the  second  program,  called  TIKIRK,  which  computes  the  intensity,  or  diffraction, 
pattern  of  the  transmitted  beam  in  the  far  field.  The  third  program,  called 
DISPLAY,  has  the  capability  of  plotting  the  temperature  and  the  intensity  results 
from  the  first  two  programs  in  a  variety  of  ways. 

The  TEMP5  heat  program  solves  the  partial  differential  heat  equation  by 
utilizing  a  numerical  integration  technique  called  the  Implicit  Alternating -Differ¬ 
ence  (I.  A.D.)  method,  in  which  all  finite  difference  analogs  of  all  derivatives  are 
correct  to  second  order.  The  method  employs  a  half  increment  shifted  net  because 
this  allows  the  use  of  general  boundary  conditions.  Being  a  fairly  comprehensive 
program,  it  allows  rather  general  specifications  of  the  initial  conditions,  as  well 
as  different  boundary  conditions  of  the  linear  type  on  each  surface.  It  is  most 
accurate  when  a  large  number  of  spatial  and  temporal  points  are  employed  in  the 
computation.  The  method  is  stable  for  all  time  increments  which  may  be  changed 
under  program  control.  Various  parameters  pertaining  to  the  beam  shape  and 
window  geometry  must  be  specified  as  input,  along  with  the  material's  thermal  and 
optical  parameters  and  its  boundary  and  initial  conditions.  The  program  computes 
the  nondimens ional  temperature  rise  at  each  grid  point  in  the  window  as  a  function 
of  nondimens  ional  time,  as  well  as  two  aberration  functions  used  by  program 
TIKIRK. 

All  of  these  results,  along  with  the  remaining  material  and  configurational 
parameters  of  the  system,  are  then  fed  into  the  TIKIRK  optical  program.  Here, 
specially  designed  integration  routines  compute  the  Hankel  transforms  and  similar 
oscillatory-type  integrals  which  evolve  from  the  vector  Kirchhoff  diffraction  theory. 

2.  Bendow,  B. ,  Jasperse,  J.R.,  and  Gianino,  P.  D.  (1972)  Optics  Commun. 

5:98.  - 

3.  Bendow,  B. ,  and  Gianino,  P.  D.  (1972)  AFCRL-72-0322,  unpublished. 

4.  Bendow,  B. ,  and  Gianino,  P.  D,  (1973)  J.  of  Electronic  Mater.  2:87. 

5.  Bendow,  B.,  and  Gianino,  P.  D.  (1973)  Appl.  Optics  12:710. 

6.  Bendow,  B. ,  and  Gianino,  P.  D.  (1973)  Appl.  Phys.  2:1. 

7.  Gianino,  P.  D. ,  and  Bendow,  B.  (1973)  Appl.  Phys.  2:71. 


The  output  gives  the  time  behavior  of  the  transmitted  laser  beam's  intensity  in 
the  far  field. 

Program  DISPLAY  is  a  general  purpose  program  which  may  be  used  to  plot 
any  function  of  two  variables  in  any  one  or  all  of  three  ways  —  contour  map, 
perspective  view  and  multiple  x-y  plot.  Since  the  temperature  in  the  window  and 
the  intensity  in  the  far  field  can  be  regarded  as  being  functions  of  two  space  vari¬ 
ables  (as  will  shortly  be  shown),  these  lend  themselves  quite  readily  to  plotting 
by  this  program.  The  third  variable,  time,  is  then  considered  to  be  a  parameter, 
as  far  as  the  two-dimensional  displays  are  concerned.  As  mentioned  above,  this 
program  has  the  very  desirable  capability  of  a  perspective  plot.  In  this  case,  a 
surface  represendng  temperature  or  intensity  as  a  function  of  the  two  space 
variables  can  be  drawn,  giving  the  appearance  of  being  a  three-dimensional  plot. 

2.  PLAN 

We  have  divided  this  report  into  two  parts.  Our  objective  in  Volume  I  is  to 
give  a  brief  introduction  to  the  problem  and  to  present  the  heat  and  diffraction 
equations  appropriate  to  the  geometric,  physical,  boundary  and  initial  conditions 
imposed  upon  the  problem.  In  this  portion,  we  also  define  all  of  the  parameters 
associated  with  the  beam  and  the  window  and  describe  the  four  common  types  of 
boundary  conditions  expected  to  occur  in  the  problem.  Finally,  we  present  in 
graphical  form  the  results  of  a  model  problem  in  which  a  typical  laser  beam  is 
incident  on  both  a  disc-  and  annular -shaped  ZnSe  window.  These  graphs  show 
various  two-  and  three-dimensional  plots  of  both  temperature  in  the  window  and 
intensity  of  the  transmitted  beam  in  space. 

Volume  II  is  a  "user's  manual."  It  describes  the  three  computer  programs  in 
detail,  explaining  how  they  function  and  how  they  are  implemented.  It  enumerates 
their  constituent  subroutines  and  subprograms,  explains  their  functions  and  gives 
Fortran  listings  in  appendices.  In  addition,  attachments  are  provided  which  give 
typical  detailed  commands  to  initiate  and  run  the  programs  in  both  the  Intercom 
and  Batch  modes  of  operation. 

3.  THE  HEAT  PROBLEM 
3.  L  The  Heat  Equation 

Since  the  incident  beam  and,  therefore,  the  temperature  distribution  in  the 
window  have  circular-cylindrical  symmetry  (there  being  no  dependence  on  the 
angular  coordinate  9),  it  will  be  sufficient  to  work  with  the  two  spatial  coordinates 


only,  viz,  r  and  z.  The  evaluation  of  the  temperature  at  any  point  in  this  window 
requires  the  solution  of  the  heat  conduction  equation,  with  the  source  term  added, 
employing  the  boundary  conditions  at  each  surface  and  taking  into  account  the 
initial  temperature  distribution  through  the  window's  interior.  The  heat  conduc¬ 
tion  equation  in  the  r,  z  cylindrical  coordinates  is  given  as: 


9T/at  =  kV,)ZT  +  Q/cpd 


(1) 


where 


T 

=  temperature 

t°C] 

t 

=  time 

[sec] 

r 

radial  coordinate 

[cm] 

z 

=  axial  coordinate 

[cm] 

d 

=  density 

[gm/cm3] 

Cp 

=  specific  heat  at  constant  pressure 

[J/gm-°C] 

K 

=  thermal  diffusivity  *  EL/c^d 

[  cm  /sec] 

K 

=  thermal  conductivity 

[W/cm-°C] 

Q 

=  volume  heating  source 

[W/cm3] 

V?. 

z  s  r-1  9/9r(r9/9r)  +  92/9z2 

[cm-2]  . 

The  window  has  an  axial  thickness  L  and  an  outer  radius  a.  When  it  is 
annular -shaped,  the  inner  radius  is  designated  by  r^  However,  rather  than  work 
with  the  coordinates  and  quantities  as  defined  above,  we  will  find  it  more  conven¬ 
ient  to  deal  with  their  dimensionless  counterparts  by  normalizing  to  the  outer 
radius  "a".  These  normalized  terms  are: 

p  =  r/a  =  radial  coordinate  (o  <  p  <  1) 

5  =  z/a  =  axial  coordinate 
Pj  =  r^/a  =  inner  radius 
p  1 2  =  (a  -  r^/a  =  radial  thickness 
£  =  L/a  =  axial  thickness 

We  choose  the  p,  £  coordinate  system  such  that  its  origin  is  at  the  center  of  the 
sample.  The  window  face  through  which  the  laser  beam  enters  is  located  at 
^  (=-L/2a).  The  laser  beam  travels  in  the  positive  ?  direction  and  exits  at  the 

face  given  by  £2  (3L/2a),  This  coordinate  system  and  the  geometrical  quantities 
defined  above  are  depicted  superimposed  on  the  window  in  Figure  1. 
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Normalized  Window  Coordinates,  "a"  is  the  outer  radius 


In  general,  the  source  term  Q  may  be  a  function  of  r  and  z  (but  not  of  time, 

except  that  it  may  be  "turned  off"  at  a  particular  time).  However,  it  is  most  con 

venient  to  assume  that  it  is  given  as  the  product  of  the  bulk  absorption  coefficient 

(3  and  the  laser  beam  power  density  I(r,  z)  at  that  point.  Typically,  I  has  dimen- 
2  ‘1 

sions  of  W/cm  and  ji  in  cm  .  The  power  density  function  is  separable  in  the 
r,  z  coordinates,  that  is,  I(r,  z)  =  I(r)I(z).  Under  the  usual  assumption  that  j3  is 
sufficiently  small  so  that  (3  times  thickness  «  1  always,  then  I(z)  «  1.  It  is  also 
assumed  that  the  incident  laser  beam  is  Gaussian-shaped  and  axisymmetric  with 
the  window.  Thus 


I(r,  z)  ~  I(r)  =  Ip  exp  (-r2/2 a2)  ,  {2' 

where  I  is  the  peak  power  density  of  the  incoming  laser  beam  and  0  is  the  beam 

radius,  ?that  is,  at  a  radial  distance  of  2a  from  beam  center,  the  power  density 

reduces  to  e~2I  .  By  dividing  a  by  the  window  radius  a,  the  normalized  beam 

radius  0  is  obtained.  Then,  Eq.  (2)  becomes  in  normalized  coordinates : 
e 

I(p)  =  Ip  exp  (-p2/2a2)  .  (3 

2 

The  relationship  between  a  ,  as  defined  here,  and  the  beam -size  parameter  a  , 

®  3 

introduced  by  Bendow  and  Gianino,  '  is  given  by: 


(4) 


The  total  power  contained  within  the  incoming  laser  beam  (P^)  is  related  to 
l(r)  in  the  following  way: 


00 

.ffwu-f 


Ip  exp  ( 


-r2/202)2wrdr  =  27ra2I 


(5) 


with  d A  being  the  elemental  surface  area.  Note  that  the  integrals  over  the  surface 
area  cover  all  space.  Solving  for  I  gives: 

Ip.Pj/2™21rP,/2  ,*\l  .  (6) 

Utilizing  Eq.  (6)  in  Eqs.  (2)  or  (3)  and  multiplying  by  (3,  we  get  the  steady-state 
volume  heating  source  term  as: 

Q(r)  =  [/3Pf/27 r02]  exp  (-r2/2 02)  , 


or. 


Q(p)  =  [PPjf /2ra2  a2]  exp  (-p2/2a2)  . 


(8) 


However,  in  most  experimental  situations  it  is  more  convenient  to  measure 
transmitted  power  (P.)  rather  than  the  incident  laser  power  (P  ).  Under  the 

t  *■  g 

assumption  that  /3L  «  1,  then  these  two  quantities  are  related  by: 

Pj./Pi  =  2n/(n2  +  1)  ,  (9) 

n  being  the  refractive  index  of  the  medium.  Solving  for  P^  : 

Pf  =  Pt(n2  +  l)/2n  .  (10) 

In  this  work  we  have  assumed  that  it  will  always  be  the  quantity  P,  that  is  known 
initially  and  the  P  will  have  to  be  computed  via  Eq.  (10).  Therefore,  the  com¬ 
puter  program  has  been  written  to  accept  P  as  input  and  to  transform  to  P^ 
according  to  the  above  formula. 

The  program  is  also  capable  of  handling  a  z-dependence  of  the  source  term  Q. 
This  capability  is  discussed  more  fully  in  Section  11.  Under  this  circumstance, 
the  assumption  of  (3L  «  1  may  not  be  valid,  necessitating  a  more  accurate  formula 
for  transmitted  power  P^  than  that  stated  in  Eq.  (9). 

3.3  The  Boundary  Condition  Equation 

As  mentioned  in  Section  3.  1,  the  window  temperature  distribution  is  also 

governed  by  the  boundary  conditions  (BC)  on  its  four  surfaces.  The  BC  Equation 

g 

for  the  ith  surface  is  of  the  form: 

KOTs/9Pi)  =  gi-hiTg.  ,  (ID 

t  l 

where  Ts.  is  the  temperature  at  the  surface,  h.  is  the  surface's  heat  transfer 

ai  n  .  »■  2 

coefficient  [W/cm  -  C],  g.  is  the  total  heat  flux  into  the  surface  [W/cm  ]  and 
d /9p.  is  the  directional  derivative  along  the  normal  pointing  outward  from  the  ith 
surface.  The  index  i  runs  from  1  to  4.  The  inner  and  outer  cylindrical  surfaces 
located  at  r1  and  "a"  are  designated  by  the  indices  1  and  2,  respectively;  whereas, 
the  plane  window  faces  located  at  z  =  -L/2  and  L/2  are  designated  by  the  indices  3 
and  4,  respectively.  There  will  be  an  equation  of  type  (11)  for  each  surface  of  the 
window,  thereby  allowing  independent  treatment  of  each  boundary. 

8.  Weil,  R.  (1970)  J.  Appl.  Phys.  41:3012;  and  Bendow,  B.,  Hordvik,  A., 

Lipson,  H. ,  and  Skolnik,  L.  (1972)  AFCRL-72-0404,  unpublished,  p.  12. 

9.  Carslaw,  H.  S. ,  and  Jaeger,  J.C.  (1959)  Conduction  of  Heat  in  Solids.  2nd 

edition,  Oxford  Press,  London,  p.  19. 


The  terms  on  the  right-hand  side  of  Eq.  (11)  refer  to  the  total  flux  of  energy 
absorbed  by  the  surface  (g^)  and  the  total  flux  emitted  (FT).  Their  difference  is 
equal  to  the  net  flux  absorbed  by  the  surface,  given  by  K(9Ts./9pj).  The  term  g^ 
in  turn,  has  two  contributors :  the  heat  influx  from  the  surrounding  environment 
whose  ambient  temperature  is  T&  and  that  from  additional  external  heat  sources 
(F.).  That  is, 


®i 


=  h.T  + 
1  a 


F. 

1 


(12) 


There  are  4  common  BC  which  typically  characterize  practically  all  of  the 
surfaces  of  interest.  They  are: 

(BC1)  Insulated  Surface:  Since  no  radiation  can  enter  or  leave  the  surface, 
both  g.  and  In  are  set  equal  to  zero. 

(BC2)  Given  Heat  Input:  hj  is  set  equal  to  zero  so  that  there  will  be  no  heat 
flow  out  of  the  surface.  All  of  the  heat  influx  from  the  external  source  originates 
from  F.,  so  that  K(9Ts./9p.)  =  gi  =  F^. 

(BC3)  No  Thermal  Constraints:  Here,  the  surface  is  experiencing  both  an 
influx  and  efflux  of  heat  radiation.  In  general,  all  of  the  terms  in  Eqs.  (11)  and 
(12)  are  present  and  are  finite,  although  in  the  typical  case  Fi  will  usually  be  zero. 
The  heat  outflow  is  represented  by  HTc-,  in  which  TSi  is  varying.  Thus,  the  net 
flux  absorbed  is  given  typically  by  h.(T  -  Tg.).  The  condition  usually  known  as 

9  I  d  L 

Newton's  Law  occurs  when  Ts.  >  T  and  F.  =  0. 

(BC4)  Fixed  Surface  Temperature:  To  hold  the  surface  temperature  Tgi  fixed, 
h.  must  be  set  equal  to  a  large  value  in  order  that  the  surface  may  be  an  efficient 
dissipator  of  any  additional  heat  reaching  it.  Furthermore,  the  temperature  of  the 
ambient  environment  in  the  immediate  vicinity  of  the  surface  must  be  the  same  as 
that  of  the  surface  (that  is,  T&  =  TSi)  so  that  there  will  be  no  net  exchange  of  energy 
between  the  two  regions.  Therefore,  if  we  set  Fi  =  0,  then  g.  =  FTa  =  h.Ts.  and 
there  is  no  net  absorbed  flux. 

The  information  on  all  of  the  parameters  pertaining  to  these  four  BC  is  tab¬ 
ulated  in  Table  1.  Note  that  all  of  the  BC  of  practical  interest  can  be  represented 
by  an  appropriate  choice  of  g.  and  In.  This  will  allow  the  resulting  computer  pro¬ 
gram  to  be  very  flexible. 

For  the  case  in  which  the  window  is  a  full  disc,  then  g1  and  h1  (which  pertain 
to  the  nonexistent  inner  cylindrical  surface)  are  set  equal  to  zero. 
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l.  THE  FAR-FIELD  INTENSITY  PATTERN 


Originally,  the  TEMP5  (heat)  program  was  written  to  solve  the  heat  and  BC 
equations  when  they  were  in  a  nondimens ional  form.  Consequently,  its  results 
(temperatures  throughout  the  window)  would  be  expressed  in  nondimens  ional  form. 
On  the  other  hand,  the  TIKIRK  (optical)  program  was  designed  to  work  with  actual 
(that  is,  dimensioned)  quantities.  To  rectify  this  discrepancy,  the  TIKIRK  pro¬ 
gram  has  been  written  to  accept  all  of  the  normalized  output  from  TEMP5  and  to 
dimens ionalize  it  accordingly  before  using  it.  Thus,  whenever  these  two  programs 
are  to  be  used  together,  only  the  normalized  data  from  the  TEMP5  program  can 
be  utilized.  We  have  already  indicated  the  transformations  between  the  spatial 
coordinates  r,  z  and  p,  5,  The  real  time  t  and  real  temperature  T  are  related  to 
their  nondimens  ional  counterparts  t  and  w,  respectively,  via: 


a 


Table  1.  Information  on  All  the  Parameters  (g.,  h.,  F.,  T  Ts.)  for  the  Boundary 
Condition  (BC)  Equation  1 


BC 

Si 

h. 

1 

F. 

1 

Ta 

TSi 

Net  Flux  Absorbed 

1 

0 

0 

0 

- 

varies 

0 

2 

F. 

0 

finite 

varies 

F. 

1 

1 

3 

h.T  +F. 

finite 

finite 

fixed 

varies 

h.(T  -Ts .) 

1  a  1 

1  a 

4 

hiTSi 

large 

0 

T  =T-. 
a  Si 

fixed 

0 

Hi 


where 


Before  the  temperature  output  of  TEMP5  —  in  the  form  of  w(p,  £,  t)  —  is  fed 
into  TIKIRK,  it  is  integrated  over  the  thickness  of  the  window  to  obtain: 


Another  function,  F2,  is  obtained  from  FI  via: 


F2(p,  r)  =  p  2  J  Flip,  t)  p  dp  .  U7) 

o 

The  dimensionless  functions  FI  and  F2  can  be  converted  to  real  time  t  by  means 
of  Eq.  (13).  Knowing  these  two  F-functions,  one  can  determine  the  time -dependent 
aberration  functions  4>P '  ®  associated  with  the  p-  and  6-  polarized  waves.  These 
<l>  -functions  contain  all  of  the  material  properties  of  the  window  and  characterize 
its  thermal  stresses  as  functions  of  time.  They  are  of  the  form: 

<J.7(p,t)  =  a  ATc  [s]  Fl(p,t)  +  4  S7  F2(p,  t)]  ,  (y  =  p,  e)  .  (18) 

4 

The  ST  are  the  thermal  lensing  parameters,  defined  as: 


sP  =  dn/dT  +  a  n3  [(1  -  v)  P12  -  v  Pi  /2  +  a  (1  +  y)(n  -  1) 

(19) 

SP  =  -S®  =  a  n3(l  +  v)(Pn  -  P12>/8 

(20) 

Sj  =  dn/dT  +  a  n3  [p^  -  2  v  P12l  / 2  +  a(l  +  j/)(n  -  1)  , 

(21) 

where  dn/dT  is  the  temperature  derivative  of  refractive  index,  a  is  the  thermal 
expansion  coefficient,  v  is  Poisson's  ratio  and  the  p.,  are  the  Pockels  elastooptic 
coefficients.  '  These  Sj  must  be  calculated  in  advance  and  read  in  as  input  since 
the  program  does  not  compute  them. 

-The  terms  containing  the  elastooptic  coefficients  in  Eqs.  (19)  through  (21)  repre¬ 
sent  the  "isotropic  approximation"  commonly  utilized  in  the  literature  for  crystal¬ 
line  windows,  but  which  is,  strictly  speaking,  valid  only  for  amorphous  solids. 
When  the  effects  of  crystallinity  are  taken  into  account,  the  forms  of  the  above 
elastooptic  terms  can  change  significantly.  For  example,  F.  Horrigan  of 
Raytheon  Corp.  and,  later,  J.  Marburger  and  M.  Flannery  of  the  University  of 
Southern  California,  have  considered  the  case  of  a  single  crystal  window  whose 
fuv.e  is  a  (111  J  plane  and  communicated  their  somewhat  different  results  tc  us. 

The  elastooptic  terms  based  on  Horrigan's  work  have  been  reported  in  B.  Bendow, 
P.  Gianino,  Y.  Tsay  and  S.  Mitra,  Appl.  Opt.  U3,  2382  (1974),  [Eq.  (33)] .  The 
correspond iufe  Wmis  based  on  the  Ivia.  nuigOi  Fdanntry  ttiialLysis  been 

reported  in  B.  Bendow  and  P.  Gianino,  Appl.  Opt.  14,  277  (1975),  [Eq.  (3)]  . 

The  results  for  the  effective  p^'s  to  be  utilized  for  polycrystalline  media  are 
given  by  M.  Flannery  and  J.  Marburger,  Proc.  of  Fifth  Conf.  on  IR  Laser  Win¬ 
dow  Materials,  C.R.  Andrews  and  C.  L.  Strecker,  Editors,  ARPA,  Arlington, 
Virginia  (1976),  p.  781. 


Since  the  sj  have  dimensions  of  (°C)  ,  then  the  $7  will  be  expressed  in 

cm. 

9  /  c 

Previously,  '  ’  we  determined  the  transmitted  intensity  at  a  field  point  in 
space  for  an  unpolarized  beam  with  an  axisymmetric  Gaussian  distribution  incir 
dent  on  the  window.  The  field  point  in  space  is  specified  by  the  cylindrical  coord¬ 
inates  X,  p1,  in  which  X  is  the  distance  from  the  window's  center  measured  along 
an  axis  orthogonal  to  the  plane  of  the  window  and  through  its  center,  while  p'  is 
the  perpendicular  radial  distance  from  this  line  to  the  field  point.  The  Gaussian 
focal  distance,  that  is,  the  axial  distance  from  the  window's  center  to  the  Gaussian 
focal  point,  is  given  by  XQ.  These  coordinates  are  shown  in  Figure  2.  Since  the 
spatial  intensity  pattern  is  circularly  symmetrical  about  this  axial  line,  there  is 
no  dependence  on  the  cylindrical  angular  coordinate. 


FIELD  POINT  (X.P* ) 


WINDOW 


LASER 

BEAM 


Figure  2.  Cylindrical  Coordinates  for  the  Field  Point 
in  Space.  The  origin  is  at  the  center  of  the  window. 
X„  is  the  Gaussian  focal  distance 
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Sometimes,  rather  than  using  the  dimensioned  coordinates  X  and  p',  it  is 
advantageous  to  employ  the  nondimens ional  generalized  coordinates  u  and  v  in  the 
calculation  of  the  intensity.  These  coordinates  are  measured  from  the  Gaussian 
focal  point  and  are  closely  related  to  the  "optical  coordinates"  of  Born  and  Wolf. 
Specifically, 

u  =  ka2(X^  -X'1)  ,  v  =  kap'/X  ,  (22) 

where  k  is  the  wave  number  2 jr/X,  and,  X  is  the  wavelength.  (At  the  Gaussian 

focus  both  u  and  v  vanish,  by  definition,  so  X  =  Xq,  p'  =  0).  Specifically,  the 

intensity  at  a  given  point  u,  v  in  the  far  field  at  any  time  t  relative  to  the  intensity 

2  ™  6 

at  the  Gaussian  focus  at  zero  time  is  given  by: 

1 

I'(u,v,t)  =  2a4(l  -  exp  (-a2))-2  |  ^  fwfx  dp|  2  +|  fwfy  dp|  2  .  (23) 

o 

The  functions  within  the  brace  in  Eq.  (23),  which  are  to  be  integrated  over  the 


window's  face,  are  defined  as  follows: 

fw(p,  u)  =  exp  {-a2p2  -  i u p 2 / 2)  , 

(24) 

fx(p,  v,  t)  =  p  JQ(pv)  exp  (ik$p)  -  fz(p,  v,  t)  , 

(25) 

fy(p,  v,  t)  =  p  JQ(pv)  exp  (ik$S)  +  f z(p,  v,  t)  , 

(26) 

fz(p.  v,  t)  =  v”1  Jx(pv)  [exp  <ik$p)  -  exp  (ik$0)]  , 

(27) 

where  JQ  and  J1  are  the  zero-  and  first-order  Bessel  functions  of  the  first  kind, 
respectively. 

Usually,  however,  it  is  desirable  to  express  the  above  intensity  in  terms  of 
the  dimensioned  coordinates  X,  p'.  In  that  case,  the  intensity  is  determined  by: 

I(X,p',  t)  =  I'(u,v,t)  ■  X2/(X2  +  p'2)  .  (28) 

It  will  be  the  objective  of  the  TIKIRK  program  to  compute  Eqs.  (23)  through  (28). 

10.  Born,  M. ,  and  Wolf,  E.  (1964)  Principles  of  Optics.  2nd  (revised)  edition, 
Macmillan  Co. ,  New  York,  p,  437. 


Before  pursuing  the  solutions  to  the  heat  and  diffraction  problems,  it  would 
serve  our  convenience  if  we  paused  to  classify  the  input  parameters  into  categor¬ 
ies  which  describe  the  window's  geometry,  its  material  composition,  surface 
characteristics  and  the  beam  properties: 

Geometrical  Parameters:  a,  L,  r^ 

Material  Parameters:  n,  /3,  K,  c^,  d,  S? 

Surface  Parameters:  gR  In 
Beam  Parameters:  a,  P^,  X 

5.  THE  PROGRAMMED  FORMS  OF  THE  HEAT  AND  BC  EQUATIONS 
5.1  Nondimensional  Form 

The  heat  and  BC  equations  can  be  cast  into  a  nondimensional  form  by  dividing 
Eq.  (1)  by  k  AT  /a2  and  Eq.  (11)  byKATc/a,  obtaining: 

9w/9t  =  Vp  ^  w  +  (2ff2  )"  1  exp  (~p2  /  2a2)  (29) 

and 

9wg  /9(Pi/a)  =  ag./KATc  -  (ah./K)wg_  ,  (30) 

in  which  use  has  been  made  of  the  normalized  distances,  plus  Eqs.  (13),  (14),  (15), 
and  (10).  Rewriting  these  equations  in  terms  of  the  notation  used  in  the  computer 
program  (that  is,  the  quantities  printed  out),  there  results: 

9U/9(TAU)  =  vfRHO)  (ZED)  U  +  A  exp  {-  (RHO)2/2(SIG)2 }  (31) 

9USi/9rii  =  G1(i)  '  H1(i)  '  Us.  (32) 

where 

V2rho)  (zed)  s  (RHO)"1  9/9(RHO)  {(RHO)  9/9(RHO)}  +  92/9(ZED)2  . 

The  integration  on  RHO  proceeds  from  RHOl  to  RH012,  while  that  on  ZED  goes 
from  -|  ZED1 1  to  |  ZEDl|. 


Upon  comparing  Eqs.  (31)  and  (32)  with  (29)  and  (30),  respectively,  one  can 
make  the  following  associations: 


TAU  h  T 
RHO  s  p 
ZED  s  £ 

RHOl  =  px 
RH012  s  p12 
ZED1  s 
A  s  1/2  <s\ 


Gl(i)  s  ag./KATc 

Hl(i)  =  ah./K  . 

rj.  =  p^a  and  refers  to  the  normalized  coordinates  p  or  £,  depending  upon  the 
surface  being  considered.  We  emphasize,  once  again,  that  the  computer  variable 
quantities  on  the  left  of  Eqs.  (33)  are  those  printed  out  by  the  program,  while 
those  quantities  on  the  right  represent  the  interpretations  that  are  to  be  given  to 
the  computer  variables.  Note  that  all  of  these  quantities  are  normalized  so  as  to 
be  dimensionless.  The  interpretation  delineated  in  Eq.  (33)  is  mandatory  when 
the  TEMP5  program  is  to  be  used  as  a  predecessor  to  the  TIKIRK  program  or 
when  temperature  plots  are  to  be  made. 

A  perusal  of  the  right  hand  sides  of  Eqs.  (33)  shows  that  one  must  know  all  of 
the  geometrical,  material,  surface  and  beam  parameters  listed  in  the  previous 
section  (except  sj  and  \)  in  order  to  specify  completely  the  computer  variables  on 
the  left  sides  of  Eqs.  (33).  For  the  case  of  no  volume  source,  the  last  term  in 
the  heat  equation  vanishes,  obviating  the  need  for  SIG  and  A.  (A  control  to  bypass 
this  term  is  available  in  the  program. ) 

5.2  A  Dimensioned  Form 

For  those  situations  in  which  one  wants  only  the  numerical  results  of  the  real 
temperature  distribution  throughout  the  window  as  a  function  of  real  space  and  time, 
but  not  the  diffracted  laser  beam  intensity,  a  different  arrangement  of  the  heat  and 
BC  equations  can  be  formulated,  resulting  in  a  more  convenient  (dimensioned) 
interpretation  of  the  computer  variables.  Here,  Eq.  (1)  is  divided  by  k  and  Eq. 

(11)  by  K,  obtaining: 


.  ‘  ■ 


3T /9t'  =  V2  z  T  +  (ATc/2ct2)  exp  <-r2/2o2) 


9Ts./9Pi  =  «i/K  ‘  <VK)  TSi  ' 

in  which  Eqs.  (15)  and  (9)  have  been  used,  and, 


t1  =  Kt  [cm  ]  . 

Again,  Eqs.  (34)  and  (35)  must  be  compared  to  the  computerized  forms, 
Eqs.  (31)  and  (32),  respectively,  resulting  in  the  following  new  associations  for 
the  printed  output: 


u 

= 

T 

[°C] 

TAU 

= 

T1 

r  2. 

[cm  ] 

RHO 

s 

r 

[cm] 

ZED 

& 

z 

[cm] 

RHOl 

= 

rl 

[cm] 

RHO  12 

= 

a  -  rx 

[cm] 

ZED1 

s 

-L/2 

[cm] 

A 

= 

ATc/2a2 

[°C/cm2 

SIG 

= 

a 

[cm] 

Gl(i) 

= 

gt/K- 

[°C/cm] 

Hl(i) 

= 

h./K 

[cm'1] 

rj.  =  p.  and  refers  to  either  r  or  z,  depending  upon  the  surface  being  considered. 
Again,  the  terms  on  the  left  of  Eqs.  (37)  are  the  computer  variables,  while  those 
on  the  right  are  their  physical  interpretations.  From  an  inspection  of  Eqs.  (37), 
it  is  obvious  that  the  temperature  as  well  as  all  distances  and  lengths  are  given  in 
their  accustomed  units;  whereas  the  time  (TAU),  volume  heating  source  (A)  and 
the  surface  parameters  (Gl(i)  and  Hl(i))  are  not. 

Once  again,  a  perusal  of  the  right  hand  sides  of  Eqs.  (37)  points  out  that  one 
must  know  all  of  the  required  parameters,  except  sT  and  X,  in  order  to  specify 
completely  the  computer  variables  on  the  left  sides  of  Eqs.  (37).  As  before,  if 
there  is  no  volume  source,  SIG  and  A  also  become  unnecessary. 
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6.  GRAPHS  FOR  TWO  MODEL  PROBLEMS 


6.1  Description  of  the  Models 

The  DISPLAY  program  has  the  capability  of  exhibiting  the  distributions  of  the 
temperature  (T)  and  the  intensity  (I)  in  a  variety  of  ways.  For  both  T  and  I  the 
program  can  provide  three  different  types  of  displays,  viz,  plots  of  contour  maps, 
perspective  views  which  give  a  three-dimensional  effect,  and  multiple  cross- 
sections  through  these  perspective  views  which  are  parallel  to  the  coordinate  axes. 
In  this  section,  we  present  various  examples  of  each  kind  of  plot  for  two  model 
problems.  These  plots  were  produced  on  the  AFGL  Calcomp  plotter  associated 
with  the  CDC6600  computer. 

The  first  model  chosen  (Model  No.  1)  consisted  of  a  disc-shaped  cylindrical 
window  of  ZnSe.  It  was  1.  3  cm  thick  and  had  a  radius  of  30  cm.  The  entrance 
face  (at  z  =  -0.  65  cm)  and  the  outer  rim  were  maintained  at  a  fixed  temperature  of 
0°C,  while  the  exit  face  (at  z  =  +0.  65  cm)  was  held  at  -100°C.  Thus,  the  three 
surfaces  satisfied  BC  No.  4  (see  Table  1).  The  remainder  of  the  window  was 
initially  at  0°C  in  an  ambient  atmosphere  also  at  0°C.  In  order  to  insure  adequate 
removal  of  heat  at  the  three  surfaces  so  that  these  temperatures  could  be  main¬ 
tained,  the  surfaces'  heat  transfer  coefficients  (h.;i=2,  3,  4)  were  set  equal  to  the 
arbitrarily  large  value  of  6  X  10A  W/cra  -  C.  From  Table  1,  the  total  surface 
flux  is  gt  =  hjTg^  Therefore,  g2  =  gg  =  0,  while  g4  =  h4Ts4  =  -6  x  1014  W/cm2. 

The  normalized  radius  of  the  laser  beam  irradiating  the  window  was  given  by 

2 

og  =  0.  289  (corresponding  to  a  =  3). 

The  second  model  (Model  No.  2)  consisted  of  an  annular -shaped  cylindrical 
window  of  ZnSe,  having  the  same  thickness  and  outer  radius  as  that  in  Model  No.  1, 
but  with  an  inner  radius  of  7.  5  cm.  The  same  laser  beam  was  assumed  here  also. 
The  inner  edge  of  this  annulus  was  fixed  at  -100°C.  Its  surface  heat  transfer 
coefficient  was  set  at  104  W/cm2  -  °C,  with  g^  =  -10®  W/cm2  (BC  No.  4).  The 
other  three  surfaces  were  left  free  to  radiate  in  accordance  with  Newton's  law, 
each  having  transfer  coefficients  of  10  2  W/cm2  -  °C  (BC  No.  3).  The  remainer 
of  the  window  was  initially  at  0°C.-  The  ambient  atmospheric  temperature  about 
the  latter  three  surfaces  was  0°C,  resulting  in  g2  =  g3  =  g4  =  0;  while  that  in  the 
vicinity  of  the  inner  surface  was  -100°C. 

6.2  Multiple  X-Y  Plots 

6.  2.  1  X-Y  TEMPERATURE  PLOTS 

For  multiple  plots  involving  the  temperature,  one  may  graph  temperature 
either  as  a  function  of  (dimensioned)  axial  distance  (z)  through  the  window  at  a 
constant  (dimensioned)  radial  position  (r),  or,  as  a  function  of  radial  distance  at 
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a  constant  axial  position.  We  arbitrarily  classify  the  first  type  as  a  "Y-plot, "  the 
second  as  an  "X-plot. " 

Examples  of  both  kinds  of  these  plots  are  shown  in  Figures  3  to  11.  Figures 
3  to  6  are  of  the  Y-type  for  the  case  of  Model  No.  1;  Figures  7  to  11  are  of  the 
X-type,  of  which  Figures  7  to  9  pertain  to  Model  No.  1,  while  Figures  10  and  11 
pertain  to  Model  No.  2.  On  all  of  these  X-type  plots,  the  (dimensioned)  radial 
distance  is  labelled  "RHO."  Also  included  in  the  caption  of  each  figure  is  the 
appropriate  plot  command  required  to  construct  the  graph.  This  plot  command 
will  be  explained  fully  in  Volume  II,  Section  10.  Figure  3  is  an  example  of  T 
versus  z  through  the  center  line  of  the  window,  that  is,  r  =  0,  for  the  9  times 
listed  on  the  right  hand  side  of  the  graph.  To  the  left  of  the  graph  are  listed 
many  of  the  TEMP5  parameters  of  interest  (see  Table  2,  Volume  II).  All  of  these 
parameters  involving  RHO's,  ZED'S,  TAU's,  Gl's  and  Hi's  are  the  normalized 
terms  enumerated  on  the  left  hand  sides  of  Eqs.  (33);  whereas,  the  quantities 
from  BETA  and  below  are  listed  with  their  actual  dimensioned  values. 

In  Figures  4  and  5,  T  versus  z  is  plotted  at  constant  radial  distances  of  15 
and  30  cm,  respectively,  for  the  same  9  times  as  in  Figure  3. 

Figures  3,  4,  and  5  are  examples  of  cases  in  which  r  was  held  fixed  for  each 
given  plot.  Figure  6  is  an  example  of  how  the  parametric  values  of  r  can  be 
varied  within  one  given  plot.  Here,  T  is  plotted  against  z  for  the  3  different  rad¬ 
ial  positions  of  0,  15  and  30  cm  at  the  fixed  time  of  5  sec. 

In  Figures  7A  and  B  the  temperature  is  plotted  against  r  in  21  different  axial 
positions  (or  planes)  for  time  values  of  1  and  5  sec,  respectively.  Two  of  these 
planes  are  singled  out  in  Figures  8  and  9,  viz,  the  one  which  is  0.  338  cm  to  the 
left  of  the  window's  center  and  the  one  0.  307  cm  to  the  right,  respectively.  These 
planes  occur  approximately  1/4  and  3/4  of  the  way  through  the  window,  respec¬ 
tively.  In  both  graphs,  T  is  plotted  against  r  for  the  9  times  shown. 

Figures  10  and  11  are  quite  similar  to  Figures  8  and  9,  except  that  Model 
No.  2  (that  is,  the  annular  window)  was  used  and  that  the  two  planes  selected 
occurred  approximately  at  the  window's  center  and  at  the  exit  face,  respectively. 
Many  of  the  TEMP5  parameters  of  interest  are  enumerated  on  the  left  side  of 
Figure  10. 

6.2.2  X-Y  INTENSITY  PLOTS 

For  multiple  plots  involving  the  intensity,  one  may  plot  intensity  either  as  a 
function  of  axial  distance  in  the  far  field  at  a  constant  radial  position,  or,  as  a 
function  of  radial  distance  at  a  constant  axial  position.  (See  Figure  2  for  a  clarifica¬ 
tion  of  axial  and  radial  coordinates  in  the  far  field).  A  control  exists  in  the  program 
so  that  the  axial  and  radial  coordinates  referred  to  here  may  be  either  the  dimen¬ 
sioned  quantities  X  and  p',  or,  the  nondimensioned  quantities  u  and  v  (see  Section  4). 
As  before,  we  classify  the  I-versus-X(or  u)  type  of  graph  as  a  "Y-plot"  and  the 
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Figure  10,  Temperature  Versus  Radial  Distance  in  the  Window  in  an  Axial  Plane  Approximately  Through  the  Window's 
Center  at  the  9  Times  Indicated.  Model  No.  2  was  used.  Plot  command  is  PLOT  (1,  1,  X,  100,  11,  NAS) 
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Figure  11.  Same  as  Figure  10,  Except  That  the  Plane  is  Located  Just  Inside  the  Window 
Plot  command  is  PLOT  (1,  1,  X,  100,  21,  NAS) 


I-versus-p'(or  v)  as  an  "X-plot.  "  However,  this  latter  X  is  only  a  code  symbol  and 
should  not  be  confused  with  the  cylindrical  coordinate  X  depicted  in  Figure  2.  The  code 
symbols  X  and  Y  are  used  only  in  the  plot  command  (refer  to  Volume  II,  Section  10). 

Examples  of  both  kinds  of  these  plots  are  shown  in  Figures  12  to  14,  in  which 
all  distances  are  measured  in  cm.  Figure  12  is  of  the  Y-type,  while  Figures  13 
and  14  are  of  the  X-type.  All  plots  pertain  to  Model  No.  1. 

In  Figure  12  the  intensity  distribution  along  the  axial  coordinate  itself  (that 
is,  zero  radial  distance)  is  depicted  for  the  9  times  indicated.  The  quantities 
cataloged  in  the  first  column  on  the  left  hand  side  of  the  graph  are  the  same 
TEMP5  parameters  as  displayed  in  Figure  3;  those  appearing  in  the  upper  portion 
of  the  second  column  are  the  relevant  TIKIRK  parameters  (see  Table  4,  Volume 
ID. 

The  intensity  along  the  radial  direction,  perpendicular  to  the  axial  coordinate 
and  at  a  fixed  distance  of  800  m  along  the  axis,  is  exhibited  in  Figure  13  for  the 
9  times  indicated.  In  Figure  14  the  same  kind  of  information  is  shown  in  the 
Gaussian  focal  plane  after  a  time  of  3  sec.  At  this  particular  time,  this  focal 
plane  is  occurring  at  a  distance  of  830  m  along  the  axis. 

6.3  Perspective  Plots 

Perspective  plots  are  intended  to  give  the  viewer  a  three-dimensional  aspect 
of  the  physical  quantity  being  observed.  Each  perspective  plot  consists  of  two 
parts:  a  calibrated  rectangular  frame  (to  indicate  the  ranges  of  the  distances  being 
covered)  with  an  arrow  (to  indicate  the  viewing  angle),  and,  the  3D  plot  itself. 

The  two  "horizontal"  coordinates  of  the  3D  plot  pertain  to  the  two  orthogonal  dis¬ 
tances  involved;  the  "vertical"  coordinate  pertains  to  the  magnitude  of  the  physical 
quantity. 

We  now  show  perspective  plots  for  both  temperature  in  the  window  and  inten¬ 
sity  in  the  far  field. 

6.  3.  1  PERSPECTIVE  TEMPERATURE  PLOTS 

The  first  part  of  the  perspective  temperature  plot  exhibits  a  calibrated  rec¬ 
tangular  frame  representing  the  transverse  rectangular  cross-section  through  the 
center  of  the  window  (that  is,  the  plane  of  the  cross-section  is  perpendicular  to  the 
window's  faces).  The  borders  of  the  frame  are  parallel  to  the  p  and  £  axes. 

Because  of  the  rotational  symmetry,  only  one  half  of  the  section  need  be  shown. 

The  horizontal  coordinate  of  the  frame  refers  to  radial  distance  along  the  window 
starting  from  its  center  (at  the  left  end  of  the  frame)  and  going  to  its  edge  (at  the 
right  end).  The  vertical  coordinate  represents  the  axial  distance  through  the 
window,  starting  from  the  entrance  face  (bottom  horizontal  line)  and  going  to  the 
exit  face  (top  horizontal  line).  The  arrow  above  the  frame  shows  the  angle  at 
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Figure  12.  Laser  Beam  Intensity  in  the  Far  Field  Versus  Axial  Distance  Along  the  Axial  Coordinate  Itself  (Zero  Rad 
ial  Distance)  at  the  9  Times  Indicated,  Using  Model  No.  1.  Plot  command  is  PLOT  (1,  1,  Y,  100,  1,  NAS) 


which  the  observer  is  to  view  the  ensuing  3D  plot  relative  to  the  rectangular 
frame.  The  angle  of  the  arrow  is  measured  counterclockwise  relative  to  the  hor¬ 
izontal  axis. 

The  second  part  of  each  perspective  plot  displays  the  actual  3D  plot  of  tem¬ 
perature  versus  radial  and  axial  distances. 

In  Figure  15A  we  see  the  rectangular  coordinate  frame  for  temperature 
appropriate  to  Model  No.  1.  Some  of  the  data  pertaining  to  the  surface  normali¬ 
zations  is  shown  in  the  legend  above  the  frame.  The  radial  distance  (RHO)  varies 
from  0  to  30  cm  as  the  axial  distance  (Z)  runs  from  -0.  65  to  +0.  65  cm.  The 
viewing  angle  is  45°.  Thus,  the  laser  beam  can  be  considered  as  travelling 
vertically  upward  through  the  frame,  \vhile  the  observer  is  on  the  exit  side  of  the 
window  looking  back  towards  the  window  at  an  angle  of  45°  with  respect  to  the 
plane  of  its  exit  face. 

The  3D  plots  themselves  are  shown  in  Figures  15B,  Q^nd  D  for  times  corre¬ 
sponding  to  10'10  (equivalently,  0),  land  8  sec,  respectively.  The  vertical 
height  of  each  3D  plot  is  proportional  to  temperature.  The  corner  of  the  plot 
closest  to  the  reader  corresponds  to  the  arrowed  corner  in  Figure  15A.  The 
narrow  horizontal  axis  represents  the  outer  edge  itself,  while  the  wider  horizontal 
axis  represents  the  exit  face  of  the  window.  As  time  increases  (Figures  15B 
through  D),  we  see  that  the  temperature  in  that  portion  of  the  window  closest  to  the 
exit  face  decreases.  This  is  in  keeping  with  the  fact  that  that  portion  of  the 
window  was  initially  at  0°C,  but  since  it  is  in  proximity  to  a  highly  radiating  sur¬ 
face  maintained  at  -100°C,  its  heat  flow  outward  is  greater  than  the  heat  input,  so 
its  temperature  will  diminish  with  time. 

A  rectangular  frame  is  always  drawn  with  each  3D  plot.  However,  since  the 
rectangular  frames  corresponding  to  the  rest  of  the  3D  plots  in  Figure  15  are 
identical  to  that  shown  in  Figure  15A,  they  are  omitted  here. 

6,  3.  2  PERSPECTIVE  INTENSITY  PLOTS 

Just  as  in  the  case  of  the  perspective  temperature  plot,  the  first  part  of  the 
perspective  intensity  plot  exhibits  a  dimensioned  rectangular  frame  giving  the 
ranges  of  the  distances  along  the  axial  and  radial  directions.  The  frame's  vertical 
coordinate  refers  to  distance  from  the  window  along  the  axial  direction;  the  hori¬ 
zontal  coordinate  to  radial  distance  measured  orthogonal  to  this  direction.  The 
arrow  above  the  frame  indicates  the  viewing  angle,  which  is  defined  exactly  as  in 
the  previous  section. 

The  second  part  of  each  perspective  plot  displays  the  actual  3D  plot  of  intensity 
versus  radial  and  axial  distances. 

In  Figure  16A  we  see  the  rectangular  coordinate  frame  for  intensity  appropri¬ 
ate  to  Model  No.  1,  together  with  its  appropriate  normalizing  information.  The 
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Figure  15.  Perspective  Temperature  Plot  for  Model  No.  1.  In  (A)  are  shown  the 
rectangular  coordinate  frame,  indicating  the  ranges  of  the  radial  and  axial  distances 
through  the  window  for  the  ensuing  3D  plot,  and,  the  arrow  indicating  the  viewing 
angle?  In  (B),  (C)  and  (D)  are  the  3D  plots  for  times  10"lu  (that  is,  0),  1  and  8 
sec,  respectively.  Plot  command  is  PERSPECTIVE 
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Figure  16.  Perspective  Intensity  Plots  for  Model  No.  1.  In  (A)  are  shown  the 
rectangular  coordinate  frame,  indicating  the  axial  and  radial  ranges  in  the  far 
field,  and,  the  arrow  indicating  the  viewing  angle.  In  (B),  (C),  (D),  (E)  and  (F 
are  the  3D  plots  for  times  10“ TO  (that  is,  0),  1,  2,  3  and  8  sec,  respectively. 
Plot  command  is  PERSPECTIVE 
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axial  distance  (X)  runs  from  500  to  1100  m  from  the  window,  as  the  radial  dis¬ 
tance  (RHO-PRIME)  varies  from  0  to  10  cm.  Again,  the  viewing  angle  is  45°. 

The  3D  plots  are  shown  in  Figures  16B  through  16F,for  times  corresponding 
to  10-10  (that  is,  0),  1,  2,  3,  and  8  sec,  respectively.  The  vertical  height  of 
each  plot  is  proportional  to  intensity.  The  corner  of  the  plot  closest  to  the  reader 
corresponds  to  the  point  X  =  1100  m,  RHO-PRIME  =  10  cm.  The  narrow  horizon¬ 
tal  axis  represents  the  axial  distance,  while  the  wider  horizontal  axis  represents 
the  radial  distance.  As  time  increases  (Figures  16B  through  1610,  we  see  the 
effects  of  thermal  lensing  (that  is,  distortion  and  defocusing)  in  the  transmitted 
beam.  In  the  case  of  ZnSe,  the  diffraction  maximum  shifts  in  towards  the  window 
(away  from  the  reader). 

The  rectangular  frames  corresponding  to  the  rest  of  the  3D  plots  in  Figure  16 
are  also  omitted  since  they  are  all  identical  to  the  one  in  Figure  16A. 

A  control  is  available  in  the  program  giving  the  programmer  the  option  of 
displaying  all  of  these  intensity  plots  in  terms  of  the  normalized  coordinates  u  and 
v,  rather  than  the  dimensioned  quantitive  X  and  p1. 

6.4  Contour  Plots 

Contour  plots  are  designed  to  exhibit  contours,  or  "profiles,  "  of  constant 
values  of  temperature  through  the  window,  or  of  intensity  in  the  far  field.  The 
same  rectangular  frames  described  in  Sections  6.  3.  1  and  6.  3,  2  are  used  and  the 
contours  are  superimposed  upon  them. 

A  set  of  temperature  contours  are  shown  in  Figures  17A  and  B  at  times  of 
2  and  8  sec,  respectively,  for  Model  No.  1.  The  symbols  designating  20  different 
temperatures  are  listed  on  the  left.  Since  the  top  horizontal  line  of  the  frame 
represents  the  exit  surface  of  the  window  (which  is  maintained  at  -100°C),  the 
lower  temperatures  occur  in  the  upper  half  of  the  frame.  Figure  17B  demonstrates 
that  as  time  increases,  the  low  temperature  contours  move  toward  the  entrance 
face  of  the  window. 

We  should  like  to  point  out  that  Figures  15D  and  17B  are  displaying  tempera¬ 
tures  for  identical  circumstances,  that  is,  the  temperature  in  the  window  for 
Model  No.  1  after  8  sec.  Moreover,  the  reader  should  keep  in  mind  that  the  cor¬ 
ner  of  the  perspective  plot  closest  to  him  in  Figure  15D  corresponds  to  the  upper 
right  hand  corner  of  the  frame  in  Figure  17B. 

Contour  plots  of  20  intensity  values  for  Model  No.  1  at  times  of  2  and  8  sec 
are  shown  in  Figures  18A  and  B,  respectively.  These  two  plots  can  be  matched  up 
with  the  perspective  intensity  plots  of  Figures  16D  and  F,  respectively.  As  above, 
the  corner  of  the  3D  plot  closest  to  the  reader  corresponds  to  the  upper  right  hand 
corner  of  the  rectangular  frame.  In  Figure  18A  we  see  that  the  diffraction 


maximum,  having  a  value  of  1.01,  occurs  approximately  860  m  from  the  window 
along  the  axial  direction.  In  Figure  18B,that  diffraction  maximum  has  diminished 
to  a  value  of  0.  876  and  is  only  about  790  m  from  the  window.  Note  in  Figure  16F 
(and  somewhat  in  Figure  18B)  how  a  subsidiary  maximum  of  almost  the  same 
intensity  value  is  occurring  at  approximately  700  m  along  the  axis,  and  how  a 
minor  peak  is  developing  approximately  3  cm  off  the  axis  at  a  distance  of  about 
950  m. 
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